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Technical  Report  Summary 

P-wave  ground  notion  spectra  at  teleselsmlc  distances  have  been 
obtained  for  several  nuclear  events.  These  include  several  events  In 
tuff  and  rhyolite,  and  one  in  granite,  at  MTS,  several  USSR  events  at 
the  Kazakh  test  region-,  and  one  event  in  the  Aleutians.  The  observed 
spectra  for  these  events  were  compared  to  theoretically  expected  spectra 
for  a  nuclear  explosion  source,  values  have  been  found  to  .scale  as 
the  l-sec  P-wave  spectral  amplitudes  at  both  NTS  and  Kazakh .  Similarly, 

M8  values  scale  as  the  20-sec  P-wave  spectral  amplitudes,  except  in  the 
case  of  one  NTS  event  (BOXCAR)  which  had  anomalously  large  surface  waves. 

At  long  periods,  events  with  similar  Mg  values  at  NTS  and  Kazakh  have 
similar  spectral  levels.  At  shorter  periods  (near  1  sec),  however,  the 
Kazakh  events  had  a  higher  spectral  level  than  events  at  MTS.  This 
probably  indicates  that  the  anelastic  attenuation  is  greater  beneath  MTS 
than  beneath  the  Kazakh  region.  However,  l-sec  spectral  levels  for 
P ILEDRIVER,  which  occurred  in  granite  at  KTS,  are  somewhat  higher  than 
those  for  events  in  tuff  or  rhyolite.  Consequently,  the  shot  medium  as  well 
the  Q  structure  appear  to  affect  the  spectr  ~e  at  short  periods. 

The  amplitudes  of  l-sec  P-waves  for  Aleutian  e’q.x^ions  depend  on  azimuth, 
apparently  because  of  focussing  and  defocussing  due  to  the  subducting 
Pacific  plate. 

Love  waves  were  used  to  study  nine  Eurasian  events  which  were  anomalous 
by  the  ;Mg  criterion  (using  Rayleigh  waves  for  Mg).  The  number  of 
anomalous  events  was  reduced  to  six  with  the  m^:Mg  criterion  if  Love  waves 


.V 


were  used  to  determine  Mg.  If  the  Love  wave  amplitude  spectrum  is  used, 
only  four  of  the  earthquakes  were  found  to  be  anomalous. 

Complex  rupture  processes  are  being  studied  by  comparing  observed 
teleseismic  P-waves  on  long-period  records  to  those  obtained  by  computing 
theoretical  seismograms  using  the  Haskell-Thomson  method.  Two  events 
were  selected  which  occurred  in  a  region  of  well  known  crust  and  upper 
mantle  structure.  The  P-wave  for  one  event  can  be  synthesized  very  well 
with  a  simple  rupture  process  (a  single  event)  whereas  the  other  event 
appears  to  consist  of  two  or  more  separate  ruptures. 

Observed  Rayleigh  waves  from  earthquakes  with  known  fault-plane 

solutions  have  been  compared  to  theoretical  Rayleigh  wave  seismograms  in 
an  attempt  to  assess  the  effects  of  propagation  path  and  near-source 
complexities  on  the  wave  forms  of  Rayleigh  waves.  For  earthquakes  which 
occur  in  the  transition  region  between  the  stable  crust  of  northern 
Asia  and  the  Himalayas,  variations  in  near  source  structure  do  not 
appear  to  have  a  great  effect  on  the  gross  features  of  Rayleigh  wave 
forms.  Minor  departures  from  theoretically  expected  wave-forms  occur 
in  almost  all  cases,  but  it  is  difficult  to  separate  the  effects  of 
near-source  complexities  from  those  of  multi-pathing  or  departures  from 
an  ldealllzed  source  geometry.  Observed  Rayleigh  wave  forms  for  an 
event  on  Novaya  Zemlya  differ  greatly  from  those  for  expected  synthetic 
seismograms.  This  could  be  due  to  complex  near-source  structure,  the 
presence  of  a  portion  of  oceanic  path  near  the  source,  or  a  source  spectrum 
which  differs  greatly  from  that  which  would  be  expected  from  presently 
accepted  models  for  nuclear  explosions. 
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INTRODUCTION 


Research  for  this  report  period  has  been  concerned  with  the  effect 

which  various  propagation  effects  and  source  properties  have  on  the 

spectra  and  wave  forms  of  body-  and  surface-waves.  The  research  is 

divided  into  four  different  but  related  reports  which  are  presented 

in  the  following  sections.  The  first  report>by  Otto  W.  Nuttli  and 

John  A.  Wagner  Vs  concerned  with  the  scaling  of  teleseismic  P-wave 

ground-motion  spectra  of  underground  nuclear  explosions.  The  second, 

■  "  - - —  — 

by  Huei-Yuin  Wen ^describes  the  use  of  Love  waves  as  an  explosion- 

t —  “  - - 

earthquake  discriminant.  The  third,  by  Jaime  Yamamoto,  describes  a 

¥ 

method  for  modelling  complex  earthquakes  and  presents  comparisons 
with  observed  data.  The  fourth  Brian  J.  Mitchell  and  Chiung-Chuan 
Cheng,  reports  work  on  assessing  the  effect  of  complex  near-source 
structure  on  the  wave  forms  of  Rayleigh  waves. 


md/^ 

-r 
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THE  SCALING  OF  TELESEISMIC  P-WAVE  GROUND-MOTION 
SPECTRA  OF  UNDERGROUND  NUCLEAR  EXPLOSIONS 

by 

Otto  W.  Nuttli  and  John  A.  Wagner 
Introduction 

Body-wave,  ,  and  surface-wave,  Mg,  magnitudes  of  underground 
nuclear  explosions  have  been  used  to  discriminate  them  from  earthquakes 
and  to  estimate  their  yield.  To  arrive  at  m^  and  Mg  values,  amplitude 
measurements  are  made  in  the  time  domain  at  periods  of  1  and  20  sec, 
respectively.  Spectral  measurements,  on  the  other  hand,  provide  amplitude 
information  over  a  broad  range  of  frequencies  rather  than  at  two  discrete 
,  .  points  and  thus  potentially  carry  more  information.  Furthermore,  the 

vf 

explosion  spectra  may  be  thought  of  as  the  limiting  case  of  earthquake 
spectra  for  earthquakes  having  a  small  source  volume  and  time  duration, 
and  thus  the  scaling  of  explosion  spectra  should  provide  a  limit  for 
earthquake  spectral  scaling,  about  which  much  remains  to  be  learned. 

This  report  is  concerned  with  an  investigation  of  the  teleseismlc 
P-wave  ground  motion  spectra  of  explosion-generated  waves.  We  intend 
to  follow  it  up  with  an  investigation  of  the  spectra  of  earthquakes, 
including  those  which  are  considered  anomalous  and  non-anomalous  by 
the  :Mg  criterion,  to  see  in  what  ways  explosion  and  earthquake  P-wave 
ground-motion  spectra  differ  from  each  other. 

Far-Field  Spectra 

Murphy  (1977)  has  presented  far-field  P-wave  displacement  spectra 
as  a  function  of  yield  for  underground  nuclear  explosions  in  saturated 


tuff,  for  what  he  calls  his  modified  source  model.  These  spectra  do 
not  take  account  of  anelastic  attenuation,  whoso  effect  is  to  make  the 

spectra  drop  off  more  rapidly  at  the  shorter  periods.  Figure  1  shows  the 
P-wave  ground-motion  spectra  at  a  distance  of  30°  based  on  Murphy's 
(1977)  modified  source  model  and  three  assumed  values  of  the  Futterman  (1962) 
attenuation  operator  T/Q  (T  is  travel  time  and  Q  is  the  quality  factor). 

For  teleseismic  P  waves  Carpenter  (1967)  found  T/Q  -  1.1  and  Helmberger 
(1973)  obtained  T/Q  “  1.0.  Plotted  in  Figure  1  are  spectra  for  T/Q  *  1.0, 
0.55  and  0,  the  latter  corresponding  to  the  curves  presented  by  Murphy  (1977) 
Mote  from  the  figure  that  the  effect  of  decreasing  Q  is  to  make  the  spectral 
curves  fall  off  more  rapidly  at  the  shorter  periods.  For  example,  at  a 
period  of  1  sec  the  spectral  amplitude  for  T/Q  ■  1.0  is  more  than  an  order 
of  magnitude  less  than  that  for  T/Q  «  0.  Because  body-wave  magnitude, 
m^,  is  obtained  from  amplitude  measurements  near  1-sec  period,  it  should 
be  very  sensitive  to  the  amount  of  absorption.  Also,  for  a  given  yield 
the  "corner  period"  increases  as  the  value  of  T/Q  increases.  For  the 
curves  T/Q  *  0.55  and  1.0  the  "corner  period"  is  taken  to  be  the  period 
at  which  the  asymptotes  for  the  short-period  and  long-period  portions  of 
the  spectral  curve  intersect. 

Murphy  (1977)  used  relative  amplitudes  for  the  far-field  P-wave 
spectra.  The  absolute  amplitudes  given  on  the  ordinate  axis  in  Figure  1 
were  determined  from  P-wave  spectra  of  NTS  events  of  known  yield  which 
were  recorded  at  a  distance  of  30s.  These  spectra  are  given  in  Figure  2. 
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T/Q  «  0 
T/Q  =  0.55 

T/Q  »  1.0 
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Pigure  i .  Theoretical  P-vave  ground -motion  spectre  at  30° 
for  Murphy’s  (1977)  modified  source  model.  The  quantity  T/Q 
is  the  Puttsrman . (i962)  attenuation  operator. 
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5  v'V  Figure  ?.  Theoretical  ?-wave  ground -notion  spectra  et  30 

tj  for  5/q  -  1.0  and  observed  soectra  from  short-oeriod  seismogran 

K  at  30°  for  3  NTS  explosions. 


9 


P-Wave  Ground-Motion  Spectra  From  NTS  Events 

Springer  and  Kinneraan  (1971)  have  given  yields  of  some  NTS  events . 

Of  them  we  have  selected  HALF  BEAK  (300  let) ,  CHARTREUSE  (70  kt)  and 
DISCUS  THROWER  (21  kt)  to  compare  their  P-wave  ground-motion  spectra  at 
a  distance  of  30*  with  those  predicted  by  Morphy's  (1977)  modified  source 
model.  The  HALF  BEAK  spectrum  is  of  the  P-vave  motion  recorded  at  SCP 
(Pennsylvania)  at  a  distance  of  30*.  The  CHARTREUSE  and  DISCUS  THROWER 
spectra  are  of  the  P-wave  ground  motion  recorded  at  COL  (Alaska)  at  a 
distance  of  33*.  We  found  that  at  their  short-period  and  the  spectra  fall 
Off  amanner  closest  to  that  of  the  theoretical  T/Q  -  1.0  carve.  In 

'  j  •  •  .  .  ,  '  •  ••  ■■  : 

Figure  2  we  ahow  the  plot  of  the  three  single-station  spectra  (dashed 
lines)  end  the  theoretical  curves  of  Murphy's  model  for  T/Q  -  l.O.  The  ob¬ 
served  spectra  were  obtained  from  the  short-period  seismograms  after  cor¬ 
rection  for  instrument  response.  The  digitizing  interval 'was  0.17  sac  and 
'the  length  of  the  P-wSve  motion  sampled  was  20  sec.  The  response  of  the 
short-period  seismograph  aystem  limits  the  reliable  values  of  the  spectra- 

.  ’  ■  .  Jr. 

'  to  tha  period  rangs  of  about  0.5  to  2  or  3  sec. 


Figure  3  contains  smoothed  P-wave  ground-motion  spectre  at  an-  apicentral 
distance  of  30*  four  #  NTS  events'  detonated  in  tuff  and  rhyolite.  Table  1 
identifies  the  numbered  events.  All  the  spectre  except  those  of  HALF  BEAK 
end  BOXCAR  Were  obtained  free  the  COL  seismograms  at  a  distance  of  33* . 

,The  short-period  part  of  the  BOXCAR  spectrum  was  obtained  from  the  NNA  (Peru) 
Short -period  seismogram  at  a  distance  of  62*,  which  was  normalised  to  s 
distance  of  30*  using  the  curve  of  Langston  and  Hslmberger  (1975).  A 
portion  of  the  BOXCAR  spectrum  between  periods  of  3  to  5  sec  was  obtained 
ftoR  the,  long-period  seismogram  of  0GD  (New  Jersey),  which  ie  at  a  distance 
Of  32*..  The  HALF  BEAT  spectrum,  as  noted  previously,  was  obtained  from 
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Figure  3.  Smoothed  P-wave  ground-motion  spectra  normalized  to 
30°  for  selected  NTS  events.  The  dashed-line  portions  of  the 
soeetra  indicate  extrapolations  based  onKurphy’s  (1977)  modi¬ 
fied  source  model.  The  events  corresponding  to  the  numbers 
are  identified  in  Table  i. 


TABLE  1 


SOURCE  PARAMETERS  OP  SELECTED  NTS  EVENTS 


Number 

1 

2 

3 

4 

5 

6 
7 


Name 

Date 

®b 

MS 

BOXCAR 

Apr. 

26, 

68 

6.2 

5.^7 

-  0.22 

(51) 

HALP  BEAK 

Jun. 

30, 

66 

6.1 

4.86 

0.18 

(48) 

ZAZA 

Sep. 

27, 

67 

5.7 

4.12 

0.28 

(14) 

BRONZE 

Jul. 

23, 

65 

5.4 

3.60 

0.15 

(14) 

DURYEA 

Apr. 

14, 

66 

5.4 

3.72 

(2) 

CHARTREUSE 

May 

06, 

66 

5.4 

3.70 

0.23 

(8) 

CUP 

Mar, 

26, 

65 

— 

3.86 

0.16 

(4) 

DISCUS  THROWER  May 

27, 

66 

5.0 

3.16 

(2) 
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the  long-period  seismogram  of  OGD  (New  Jersey) ,  which  is  at  a  distance  of 

32°.  The  HALF  BEAK  spectrum,  as  noted  previously,  was  obtained  from  the 

SCP  short-period  seismogram.  All  the  spectra  of  Figure  3  were  smoothed, 

using  as  a  guide  the  Murphy  curves  for  T/Q  *=  1.0.  The  dashed  portions  of 

the  spectra  indicate  extrapolations  based  on  Murphy's  curves. 

Values  of  the  body-wave  magnitude,  m^,  and  the  surface-wave  magnitude. 

Mg,  of  the  8  events  of  Figure  3  are  given  in  Table  1.  The  m^  values  were 

taken  from  the  publications  of  the  International  Seismological  Centre. 

The  M  values  were  determined  by  us,  using  the  amplitudes  of  the  vertical- 
0 

component  Rayleigh-wave  motion  at  20-sec  period.  The  number  in  parentheses 

following  the  M  value  Indicates  the  number  of  stations  used  in  the  M 
S  b 

determination.  The  number  following  the  -  sign  is  the  standard  deviation 
of  the  average  Mg  value. 

We  can  observe  from  Table  1  and  Figure  3  that  all  the  spectral  levels 
increase  as  increases,  i.e.  the  larger  the  value,  the  larger  the 
spectral  amplitude.  In  a  more  quantitative  sense,  we  can  compare  changes 
in  m^  with  changes  in  the  logarithm  of  the  spectral  level  at  1-sec  period. 

The  results  are  given  in  Table  2,  in  which  all  events  are  compared  to 
BOXCAR  (event  no.  1).  In  all  cases  except  CHARTREUSE,  the  difference  between 
Am^  andALog1^A(l  sec)  is  only  0.1  unit.  For  CHARTREUSE  it  is  0.2  units. 

This  is  truly  remarkable,  for  the  uncertainties  in  both  the  smoothed  spectral 
amplitudes  and  the  m^  values  are  at  least  0.2  units  each.  It  indicates 
that  in  spite  of  the  fact  that  the  P-wave  displacement  spectrum  is  falling 
off  rapidly  at  periods  near  1  sec,  the  m^  value  is  a  good  estimate  of  the 
spectral  amplitude  at  1-sec  period. 


Number 

1  1 

i  ■ 


3 

e  4 

fj  5 

i:  6 
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TABLE  2 

DIFFERENCES  IN  mb  AND  log1Q  A(1  sec)  FOR 
SELECTED  EVENTS,  COMPARED  TO  BOXCAR 


Name  Alog1Q  A(1  sec) 


EOXCAR 

— 

— 

HALF  BEAK 

0.1 

0.2 

ZAZA 

0.5 

0.6 

BRONZE 

0.8 

0.9 

DURYEA 

0,8 

0.9 

CHARTREUSE 

0.8 

1.0 

CUP 

— 

1.0 

DISCUS  THROWER 

1.2 

1.3 

14 


We  can  also  see  how  the  20-sec  spectral  level  scales  with  Mg.  The 
results  are  given  in  Table  3,  in  which  all  events  are  compared  to 
BOXCAR  (event  no.  1).  The  results  are  not  so  good  as  Table  2.  First 
of  all,  the  Mg  value  of  BOXCAR  must  have  been  anomalously  large,  by  about 
0.6  to  0.7  units.  This  could  result  either  from  release  of  accumulated 
tectonic  stress  by  the  explosion,  or  by  development  of  a  crack  system  with 
a  preferred  direction  of  orientation  of  the  cracks.  If  the  Mg  value  of 
BOXCAR  is  arbitrarily  taken  as  0.6  units  less  than  the  observed  value, 
then  the  differences  between AMg  and  Alog^gA(20  sec)  in  Table  3  are  no 
greater  than  0.3  units  for  the  events  studied.  Or,  to  put  it  in  another 
way,  if  HALF  BEAK  is  taken  as  the  standard  for  comparison,  then  all  the 
events  except  BOXCAR  in  Table  3  show  a  difference  between AMg  and  Al°g^gA(20  se 
of  0.3  units  or  less.  As  the  phenomenon  of  anomalously  large  surface-wave 
generation  by  release  of  tectonic  stress  or  by  extensive  cracking  is  to  be 
expected  only  for  very  large  events  (of  the  size  of  BOXCAR) ,  we  can  conclude 
for  the  somewhat  smaller  events,  of  more  interest  from  a  detection  and 
discrimination  standpoint,  that  the  20-sec  period  spectral  level  of  the  P-wave 
ground  motion  will  scale  as  Hc- 


P-Wave  Ground-Motion  Spectra  From  MILROW 
(Aleutian  Islands)  and  USSR  Explosions 


Spectra  of  the  P-wave  ground  motion  were  determined  from  single- 
station  observations  for  6  USSR  explosions  and  for  MILROW,  which  took 
place  in  the  Aleutian  Islands.  For  one  USSR  explosion  in  Novaya  Zemlya 
the  spectrum  was  determined  also  from  a  long-period  seismogram.  In 
addition  to  the  Novaya  Zemlya  event  there  were  4  in  Kazakhstan  and  one  in 


western  Russia 
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TABLE  3 


DIFFERENCES  IN 'Mg  AND  log1Q  A(20  sec)  FOR  SELECTED 
NTS  EVENTS,  COMPARED  TO  BOXCAR 


Number 

Name 

AMS 

*log10  A(20  sec) 

1 

BOXCAR 

— 

— 

2 

HALF  BEAK 

0.6 

0.2 

3 

ZAZA 

1.35 

0.8 

4 

BRONZE 

1.9 

1.1 

5 

DURYEA 

1.75 

1.1 

6 

CHARTREUSE 

1.8 

1.1 

7 

CUP 

1.6 

1.1 

8 

DISCUS  THROWER 

2.3 

1.4 

Figure  4  shows  the  smoothed  P-wave  ground-motion  spectra  of  the 
7  events ,  normalized  to  an  epicentral  distance  of  30° .  The  dashed  parts 
of  the  spectra  indicate  extrapolations  based  on  Murphy's  (1977) 
theoretical  curves,  after  they  were  modified  to  take  account  of  absorption. 
Table  4  identifies  the  events,  which  are  labeled  A  through  G.  The 
values  were  determined  by  the  International  Seismological  Centre,  and  the 
Mg  values  were  determined  by  us  from  the  amplitudes  of  20-sec  period 
Rayleigh  waves. 

From  an  examination  of  Figure  4  we  can  see  that  the  spectral 
amplitude  levels  increase  as  the  values  increase.  A  more  quantitative 
comparison  is  given  in  Table  5,  in  which  are  listed  the  m^  and  log^(l  sec) 
differences  with  respect  to  M1LR0W.  From  the  table  it  appears  that  the 
logjgA(l  sec)  value  for  MILROW  is  about  0.6  units  too  large.  If  it  were 
reduced  to  the  1-sec  amplitude  of  event  B  all  the  remaining  events  would 
scale  to  within  0.3  units.  Looking  at  Figure  4,  the  higher- than-expec ted 
spectral  level  for  MILROW  at  1-sec  period  results  from  an  apparently  lower 
average  T/Q  value  for  it  than  for  the  Novaya  Zemlya  event  B.  An  ^  value 
based  on  the  single  station  BHP  (Canal  Zone)  used  to  obtain  the  MILROW 
spectrum  also  was  found  to  be  about  0.6  unit  higher  than  the  average 
value  of  6.5.  The  effect  of  the  descending  Pacific  lithospheric  slab 
under  the  Aleutians  is  to  cause  different  apparent  absorption  for  rays 
leaving  in  different  azimuths,  and  thus  it  should  produce  different 
shaped  spectra  in  different  azimuths.  If  we  had  averaged  the  spectra 
for  the  P-wave  motion  of  MILROW  at  a  number  of  stations  we  probably  would 
have  arrived  at  a  spectrum  similar  to  that  for  the  Novaya  Zemlya  event  B. 
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Figure  4.  Smoothed  P-wave  ground -mot ion  spectra  normalized  to 
30*  for  MILHQW  and  selected  USSR  events.  The  dashed-line 
portions  of  the  spectra  indicate  extrapolations  based  on 
Murphy1 8  (1977)  modified  source  model."  The  events  corresponding 
to  the  letters  are  identified  in  Table  4. 
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TABLE  4 

SOURCE  PARAMETERS  OP  KILROW  AND  SELECTED 
USSR  EVENTS 


Letter 

A 

E 


Location 

Date 

mb 

Kg 

MILROW  -  Aleutians 

Oct.  02, 

69 

6.5 

5.09 

±  0.22 

(40) 

Novaya  Zemlya 

Sep.  27, 

71 

6.5 

5.05 

0.31 

(14) 

Kazakh 

Dec.  22, 

71 

6.0 

4.17 

0.14 

(8) 

Kazakh 

Oct.  21, 

71 

5.5 

3.43 

0.32 

(5) 

Kazakh 

Nov.  29, 

71 

5.4 

3.83 

0.51 

(5) 

Kazakh 

Oct.  09, 

71 

5.3 

3.23 

(1) 

Western  Russia 

Ocx.  04, 

71 

4.6 

<2.8 

Table  6  shows  that  the  differences  of  Mg  and  log^QA(20  sec)  with 
respect  to  MILROW  are  all  within  0.2  units  except  for  event  E,  whose 
Mg  value  had  a  very  high  standard  deviation  (0.51).  Thus  probably  the 
Mg  value  of  event  E  is  too  large.  With  that  single  exception  the  Mg  values 
scale  as  the  20-sec  period  spectral  P-wave  amplitudes.  This  indicates 
that  none  of  the  7  events  had  anomalously  large  surface  wave  generation, 
as  was  noted  for  BOXCAR. 

Comparison  of  USSR  and  NTS  Spectra 

Figure  5  compares  the  P-wave  spectra  of  some  NTS  events  (solid  lines) 
and  USSR  events  (dashed  lines),  all  normalized  to  an  epicentral  distance 
of  30*.  The  figure  shows  that  events  at  NTS  and  in  Eurasia  of  the  same 
Mg  have  similar  20-sec  P-wave  spectral  levels,  and  furthermore  that  at 
the  shorter  periods  the  USSR  spectra  indicate  somewhat  lower  absorption 
than  the  NTS  spectra.  This  latter  observation  has  two  important 
consequences:  1)  a  USSR  explosion  of  a  given  yield  will  have  a  higher 
m^  value  than  an  NTS  explosion  of  the  same  yield  (the  ra^  difference  will 
be  about  0.3  units  for  an  event  about  the  size  of  CHARTREUSE),  2)  there 
will  be  some  regional  dependence  of  the  :Mg  relation  for  explosions, 
because  for  a  given  Mg  the  m^  value  of  the  explosion  will  depend  on  the 
P-wave  absorption  in  the  source  region.  In  terras  of  T/Q,  the  average 
T/Q  for  NTS  events  is  about  1.0,  whereas  for  Kazakh  events,  it  is  about 
0.8. 

Summary  and  Conclusions 

Murphy's  (1977)  spectra  for  a  nuclear  explosion  source  model 
were  compared  with  P-wave  ground-motion  spectra  at  30*  epicentral  distance 
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DIFFERENCES  IN  Kg  AND  log1Q  A(20  sec)  FOR 
SELECTED  USSR  EVENTS,  COMPARED  TO  NILROW 


Letter 

Location 

*MS 

&log1Q  A( 20  sec) 

A 

Aleutians  (NILROW) 

— 

— 

B 

Novaya  Zemlya 

0.0 

0.1 

C 

Kazakh 

0.9 

1.1 

D 

Kazakh 

1.7 

1.5 

E 

Kazakh 

1.3 

1.7 

F 

Kazakh 

1.9 

1.8 

G 

Western  Russia 

p-2.3 

2.5 

DISPLACEMENT 
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Figure  5.  Smoothed  P-wave  ground-motion  spectra  normalized  to 
30®  for  selected  NTS  events  (solid  linos)  and  selected  USSR 
events  (dashed  lines).  Note  chat  the  USSR  explosion  spectra 
cross  the  NTS  explosion  spectra  at  the  shorter  periods,  indicat 
ing  smaller  absorption  in  the  USSR  (Kazakh)  test  sixe  region 
than  at  NTS. 
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from  NTS,  USSR  and  an  Aleutian  explosion.  For  the  NTS  and  Kazakh  test 
sites  the  log^^Ad  see),  where  A(1  sec)  is  the  amplitude  of  the  P-wave 
spectrum  at  l-sec  period ,  scaled  directly  as  ^ .  Similarly,  the 
logj^A(20  see)  scaled  directly  as  Mg,  except  for  BOXCAR  which  had 
anomalously  large  surface  waves. 

NTS  and  Kazakh  spectra  of  events  with  the  same  Mg  values  had  the 
same  long-period  P-wave  amplitude  levels.  However,  at  the  shorter  periods 
near  1  sec  the  Kazakh  spectra  had  higher  amplitudes  than  the  NTS  spectra. 
Indicating  lower  absorption  in  the  Kazakh  source  region  (average  T/Q  of 


about  0.8  for  Kazakh  and  1.0  for 


.  The  difference  amounts  to  about 


0.3  m^  units  for  an  event  the  size  of  CHARTREUSE.  Thus  a  Kazakh  explosion 
of  the  same  yield  as  an  NTS  one  will  have  a  higher  m^  value,  but  the  same 
Mg  value,  this  assumes  that  there  is  no  anomalously  large  surface  wave 
generation  due  to  the  release  of  regional  stress  or  to  extensive  cracking 
in  preferred  directions  in  the  source  region. 

The  1-see  period  amplitudes  of  P  waves  from  Aleutian  explosions  show 
s  dependence  on  azimuth,  which  apparently  is  related  to  the  subducting 
Pacific  plate.  This  is  true  for  amplitude  measurements  both  in  the 
frequency  and  the  time  domain.  The  Mg  values  do  not  show  such  an  azimuthal 
dependence ,  for  at  periods  of  20  sec  the  effect  of  the  plate  geometry 
apparently  is  insignificant . 

Finally,  the  m^:Mg  discriminant  will  show  soma  dependence  on  the 
region  of  the  explosion,  because  the  value  is  affected  by  absorption  in 
the  source  region. 


Dr.  Robert  Masse  of  AFTAC,  after  reading  a  draft  of  this  report, 
suggested  that  wa  look  at  an  NTS  shot  in  granlts.  We  selected 
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Figure  6.  Smoothed  P-wave  spectra  normalized  to  30°.  The 
solid-line  curves  are  for  the  spectra  of  PILEDRIVFR  (57  let) 
as  recorded  at  COL  and  LPS.  The  long  dashed-line  ctirves  are 
Murphy’s  (1977)  curves  for  T/Q  =1,0  and  shots  of  100  kt  and 
10  kt  yield  in  tuff/rhyolite.  The  short  dashed-line  curve  is 
for  Kazakh  event  J)  (see  T-hle  4).  The  dash-dot  curve  is  for 
NTS  events  4  and  5  in  tufi/rhyolite  (see  Table  1),  Note  that 
the  PILEDaIVFR  curves  indicate  less  apparent  absorption,  com¬ 
pared  to  tne  tuff /"hyalite  curves,  at  the  shorter  oeriods. 

The  differences  in  the  spectral  levels  of  the  ^ILPDRIVFR  ground 
motion  at  COL  and  LPS  suggest  that  absorption  in  the  crust  and 
upper  nan tie  near  the  receiver  can  also  affect  the  spectral 
levels. 
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1-0 VE  WAVES  AS  AN  EXPLOSION-EARTHQUAKE  DISCRIMINANT 

by 

Huei-Yuin  Wen 

Introduction 

The  difference  in  the  surface-wave  versus  body-wave  magnitude 
relationship  between  earthquakes  and  underground  explosions  was  first 
investigated  by  Press  et  al  (1963)  and  Brune  et  al  (1963) .  Now 
seismologists  generally  agree  that  earthquakes  and  explosions  can  be 
discriminated  above  a  certain  magnitude  level  by  using  the  ( :Mg 
criterion.  For  most  studies  (Liebermann  &  Pomeroy,  1969;  Evernden,  1969; 
Basham  et  al.  1970;  Marshall,  1970),  the  body-wave  magnitude  was 
determined  from  P  waves  of  period  about  1  sec  and  the  surface-wave 
magnitude  Mg  from  the  vertical-component  of  Rayleigh  waves  of  period, 
about  20  sec.  Evernden  (1975)  found  that  the  Mg (20)  for  explosions  vas 
1.0  to  1.25  units  lower  than  that  for  shallow-focus  earthquakes  at  the 
same  m^  (m^>  4) .  Some  anomalous  earthquakes  (earthquakes  which  have 
ra^iMg  values  nearly  those  of  explosions)  occur  near  continental  margins 
(Evernden,  1975)  but  most  occur  in  the  interior  of  the  continent  (Lander, 
1972;  Nuttli  and  Kim,  1975).  For  such  events,  the  m^Mg  criterion  of 
discrimination  appears  to  fail. 

Fundamental-mode  Love  waves  have  no  node  in  the  displacement  versus 
depth  curve.  The  amplitude  of  Rayleigh  waves  can  have  a  minimum  at 
20-sec  period  for  some  earthquakes  because  of  the  particular  depth  and 
source  mechanism,  resulting  in  anomalously  low  Mg.  In  that  case  the  Love- 
wave  magnitude  can  be  used  to  reclassify  these  earthquakes  as  non-anomalous 


ones.  Evernden  (1975)  found  that  for  underground  explosions  the  Love-wave 
amplitudes  are  significantly  less  than  the  Rayleigh-wave  amplitudes. 


whereas  for  earthquakes  they  are  of  about  the  same  value.  Therefore,  a 

L 


criterion  based  on  Love  waves  (Mg)  should  be  more  powerful  than  one 
based  on  Rayleigh  waves  In  addition,  use  of  both  an<* 

M^jM^2  criteria  should  greatly  improve  identification  capabilities. 

A  recent  study  (Nuttli  et  al.  1975)  has  shown  that  some  aftershocks 
or  foreshocks  of  non-anomalous  mainshocks  in  the  continental  interior  of 
Eurasia  possessed  anomalous  m^M^2  values.  In  this  paper,  nine  Eurasian 
events  with  anomalously  low  Rayleigh-wave  magnitudes  (from  Nuttli  et  al, 

1975)  and  fourteen  northwestern  Kashmir  mainshock-aftershock  sequence  earth¬ 
quakes  which  occurred  during  the  year  1972,  have  been  studied  to  testify 
to  the  usefulness  of  Love  waves  as  a  further  discriminant  between  earthquakes 
and  underground  explosions.  The  data  are  obtained  from  70-mm  film  copies 
of  the  VLPE  stations  and  a  selected  number  of  WWSSN  stations. 


Love-Wave  Magnitudes  (Mg) 

In  order  to  determine  M^,  one  first  must  separate  the  Love  waves 
from  the  Rayleigh-wave  motion  on  the  horizontal-component  seismograms. 

The  procedure  for  doing  this  is  as  follows. 

First,  identify  the  maximum  amplitude  wave  motion  with  period  about 
20-3  seconds  as  recorded  on  the  long-period,  horizontal-component 
seismograms  (N-S  and  E-W  components).  Second,  digitize  two  minutes  or 
more  of  the  maximum  wave  motion.  The  sampling  interval  and  the  digitization 
times  must  be  the  same  for  both  N-S  and  E-W  components.  Third,  resolve 
the  digitized  data  into  transverse  component  UT  (duo  to  Love-wave  motion) 
and  radial  component  UR  (due  to  Rayleigh-wavo  motion) . 
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After  traces  of  the  Love-wave  motion  are  obtained,  can  be  calculated 
by  the  use  of  the  1ASPE1  (INternational  Association  of  Seismology  and 
Physics  of  the  Earth’s  Interior)  formula  (B&th,  1969), 

Mj:  -  3.30  +  1.66Log1QA+  Log10-|-  (25°  *  A  <  140°)  (1) 

and  the  Nuttli-Kim  (1975)  formula, 

-  4.16  +  1.07Log1()A  +  Log10-|-  <10°  £  A  ±25°)  (2) 

where  A  (in  microns)  is  one-half  the  maximum  peak-to-peak  amplitude  of 

the  Love  waves,  T  is  the  corresponding  period  of  20  -  3  seconds  and  A 

(in  degrees)  is  the  epicentral  distance. 

The  values  for  the  nine  anomalous  earthquakes  are  given  in 

Table  1.  These  anomalous  Eurasian  events,  which  occurred  in  the  first 

half  year  of  1972,  are  not  restricted  to  any  single  geographic  area. 

Table  1  indicates  that  they  took  place  in  the  area  of  10*  -  50*N  and 

14*  -  95*E.  Events  No.  45  and  145  are  oceanic  earthquakes,  whereas  the 

others  are  continental  events.  The  focal  depths  determined  by  NEIC 

RZ 

were  all  shallow.  The  probable  factors  Influencing  the  relationship 

have  been  attributed  to  focal  depth,  focal  mechanism,  uncertainties  in 
the  crustal  thickness,  travel  path,  source-time  function  and  the  finiteness 
of  the  source  (Tsai  and  Aki,  1971;  Liebermann  and  Pomeroy,  1970;  Evernden, 
1975) . 

For  a  normal  earthquake,  Love  waves  are  at  least  as  high  in  amplitude 
as  Rayleigh  waves,  if  not  larger,  if  an  adequate  number  of  azimuths  is 
sampled.  Evernden  (1969)  found  that  some  anomalous  earthquakes  occurred 
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RZ  X 

in  Eurasia  having  larger  Love  than  Rayleigh  waves.  By  the  M  :MT 

S  0 

criterion,  the  number  of  anomalous  earthquakes  is  reduced  to  six  (events 
which  lie  above  either  the  long  or  short  dashed  lines  in  Figure  1) . 

The  northwestern  Kashmir  mainshock-aftershock  sequence  consisted  of 
14  earthquakes  with  values  in  the  range  of  4  to  6  that  occurred  during 
one  month  of  1972  and  which  had  small  variations  in  hypocentral  coordinates 
and  focal  depths  (estimated  by  NEIC) .  An  advantage  of  using  aftershocks 
is  that  it  removes  the  propagation  effect,  if  the  earthquakes  are  compared 
to  one  another  and  if  the  data  of  the  same  stations  are  used  for  all  the 
earthquakes . 

Table  2  gives  magnitudes  and  hypocentral  coordinates  of  the  Kashmir 

earthquakes.  All  earthquakes  of  the  sequence  except  for  the  mainshock  had 

0.59<mb  -  <1.26.  In  general  if  the  difference  exceeds  unity  the 

earthquakes  are  in  or  near  the  explosion  population.  These  observed 

differences  between  and  M^2  therefore  indicate  that  some  of  the  aftershocks 

have  explosion-like  ra^sM^2  values. 

Figures  2  and  3  are  the  plots  of  m^M^2  and  m^ 

RZ 

the  earthquakes.  On  the  m^:Mg  discriminant  (Figure  2),  the  mainshock 

which  occurred  on  September  3,  1972,  of  -  5.97  and  M^2  -  6.05,  appears 

to  fall  on  the  earthquake  side  of  the  curve,  whereas  the  majority  of 

aftershocks  appear  on  the  fringe  of  the  earthquake  m^rM^2  population.  The 
RZ 

anomalous  m^Mg  events,  namely  K3,  K4,  K7,  and  K12,  switch  to  the  earth¬ 
quake  part  in  the  plots,  as  seen  in  Figure  3.  All  the  events  have 

Love  waves  at  least  as  large  in  magnitude  as  the  Rayleigh  waves.  Except 
in  a  region  displaying  very  high  tectonic  stress  release,  the  ratio  of  Love- 

wave  to  Rayleigh-wave  magnitude  of  earthquakes  is  not  the  same  as  explosions, 

RZ  X 

even  multi-shot  ones  (Evernden,  1976).  Therefore,  when  both  ra^Mg  and  m^:Mg 
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values  for  all 


criteria  are  employed,  there  is  compensation  for  any  effects  of  the 
Rayleigh-wave  radiation  pattern,  and  we  can  conclude  that  the  14  Kashmir 
earthquakes  are  all  non-anomalous  earthquakes. 

Focal  Mechanisms  of  Kashmir  Earthquakes 

The  fault-plane  solutions  for  the  events  of  the  Kashmir  mainshock- 
aftershock  sequence  are  given  in  Figure  4.  The  P-wave  polarities  were  de¬ 
termined  from  the  long-period,  vertical-component  seismograms.  For  some 
stations,  the  first  arrival  of  the  P-wave  on  the  short-period,  vertical- 
component  seismograma  waa  used  when  the  long-period,  first  P-wave  motion 
was  not  available.  The  octagons  and  triangles  stand  for  compressional 
and  dilatational  P-wave  motion,  respectively. 

The  focal  mechanism  solutions  for  events  K1  and  K7  (Figures  4(a)  and 
4(g))  are  well  determined  and  are  similar  to  each  other.  They  both 
are  characteristic  of  thrust  faulting,  with  nodal  planes  striking  nearly 
in  the  north  direction.  For  the  remaining  events  a  unique  fault-plane 
solution  could  not  be  obtained  from  the  P-wave  data  alone.  But  for  these 
remaining  earthquakes  the  available  data  are  compatible  with  the  same  focal 
mechanism  for  most  of  the  Kashmir  earthquakes.  The  exceptions  are  event 
K8,  which  shows  normal  faulting,  and  event  Kll  which  has  only  two  readings 
of  P-wave  polarities  from  available  seismic  stations. 

Love-Wave  Amplitude  Spectra  (Ag) 

Because  all  the  events  are  intermediate  to  small  magnitude  events 
and  because  some  of  the  records  had  microseisms  too  high  to  detect  the 
surface  waves,  both  the  Love-wave  magnitudes  and  the  Love-wave  amplitude 


amplitude  spectra  have  been  determined  with  great  effort  and  with 
difficulty. 

The  spectral  amplitude  of  fundamental-mode  Love  waves  in  the  period 
range  10  to  50  seconds  were  computed  by  digitizing,  rotating  coordinates 
and  Fourier  analyzing  the  resultant  transverse-component  digital  data. 

The  spectra  were  corrected  for  instrument  response  and  for  linear  trend. 

They  were  then  normalized  to  a  reference  distance  of  9.0s  (1000  km)  from 
the  epicenter  by  removing  the  effects  of  geometrical  spreading,  using  a 
filter  program  (Herrmann,  1973) .  Figures  5  and  6  give  sample  spectra  for 
Eurasian  events  whose  records  could  be  digitized,  to  see  how  they  scale  over 
the  period  range  considered.  The  spectral  shape  should  also  be  compared 
for  individual  events. 

Tsai  and  Aki  (1971)  pointed  out  that  the  differences  in  amplitude 
spectral  shape  of  Love  waves  between  explosions  and  small  earthquakes  do 
not  significantly  depend  on  the  source  mechanism  and  focal  depth,  nor 
crustal  thickness.  This  implies  that  the  differences  in  the  spectral 
shape  of  Love  waves  might  be  due  to  multi-pathing,  lateral  variation  in 
velocity,  attenuation  and  radiation  pattern.  By  inspection  of  Figures 
5  and  6,  some  systematic  patterns  emerge. 

From  a  single  station  (K0N)  observation  for  event  No.  37,  the  spectral 
shape  resembles  that  of  a  normal  shallow  earthquake  (Figures  6  and  7  of 
Tsai  and  Aki,  1971);  the  Love-wave  spectrum  peaks  at  a  period  of  about 
20  sec  and  decreases  at  the  longer  periods.  Figure  5(b)  shows  Love-wave 
spectral  shapes  somewhat  unlike  that  given  for  event  No.  37.  The  spectral 
level  increases  from  20  to  40  sec,  indicating  a  somewhat  deeper  focus. 

Over  the  period  range  of  10  to  20  sec,  the  spectra  have  relatively  steep 
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short -period  decay  rates,  because  the  short-period  waves  are  absorbed 
by  the  water  layer. 

Concerning  the  remaining  events,  the  data  indicate  that  events 
No.  73,  145,  and  146  all  generate  Love  waves  with  spectral  content 
similar  to  event  No.  45,  whereas  event  No.  135  and  almost  all  the  Kashmir 
earthquakes  (Figures  5(d)  and  6)  have  shallow  focus  characteristics.  If 
the  focal  depths  of  those  events  as  determined  by  NEIC  are  correct,  it  is 
inferred  that  differences  in  the  prevailing  periods  of  surface  waves,  as 
Indicated  by  the  ratio  of  the  amplitude  spectral  density  of  20/40-sec  Love 
waves,  will  indicate  whether  these  shallow  earthquakes  are  anomalous  or 
non-anomalous  if  the  ratio  is  less  than  or  greater  than  one,  respectively. 
Table  3  gives  the  relations  between  magnitudes  and  amplitude  spectra  (in 

units  of  era-sec)  of  Love  waves.  If  the  ratio  of  Ag(20)/Ag(40)  is  assumed 
to  be  independent  of  magnitude  and  source  region,  it  can  provide  a  discrimi¬ 
nation  capability.  That  is,  for  shallow  non-anomalous  earthquakes  the 
L 

ratio  Ag 

depth  and  source  mechanism  of  the  events  in  the  Kashmir  aftershock  sequence 
are  not  different,  and  because  the  propagation  paths  are  essentially  the 
same,  the  question  as  to  how  the  discrepancy  in  the  spectral  shape  of  Love 
waves  recorded  at  station  SHI  arises  is  not  yet  answered. 

:Ag  Relationship 

From  the  concept  of  surface-wave  magnitude  (M^)  end  the  amplitude 
spectrum  (Ag) ,  Mg  should  have  a  linear  relation  with  Log^Ag  at  a  period 
of  20  sec.  In  addition,  there  is  an  equivalence  of  logarithmic  amplitude 
differences  in  the  time  domain  (^3Mg)  with  those  in  the  frequency  domain 
(dLog^Ajp.  Taking  the  data  for  event  No.  37  at  KON  as  the  reference  value 


(20)/Ag(40)  should  be  larger  than  one.  By  knowing  that  the  focal 


for  those  nine  Eurasian  events  and  SHI  for  event  K2  for  Kashmir  events. 
Table  3  shows  (4th  column  and  7th  column)  that  AM^  is  almost  consistent 
with  ALog10A£<20). 

A  relation  between  Love-wave  magnitudes  (M^)  and  amplitude  spectra 
(Ag)  at  a  period  of  20  sec  and  an  epicentral  distance  of  9*  has  been 
determined  by  a  linear  least-square  fit  to  the  data  of  Table  3.  The 
regression  line  is 


-  (0.902  ±  0.113)Log1(JAg  +  (6.688  -  0.318) 


which  is  shown  in  Figure  7,  for  93Z  confidence. 


Discussion  and  Conclusions 


Love-wave  magnitudes  and  amplitude  spectra  were  investigated  for 
Eurasian  earthquakes.  A  spectral  amplitude  determination  is  superior  to 
a  magnitude  determination  of  the  size  of  an  earthquake  because  it  covers 
a  range  of  periods.  From  the  relation  between  m^  and  and  the  relation 
between  and  Ag ,  a  hypothesis  is  proposed  that  a  more  effective  discriminant 
between  earthquakes  and  underground  explosions  is  one  based  on  a  comparison 
of  body-wave  magnitudes  (m^)  and  Love-wave  amplitude  spectra  (Ag) .  By 
combinations  of  Figures  1,  3,  and  7,  the  amplitude  spectra  data  are  plotted 
in  Figure  8  as  a  function  of  body-wave  magnitude  for  all  the  Eurasian 
events  studied  in  the  paper.  This  figure  shows  that  only  four  events 
(events  No,  45,  73,  145,  146)  are  classified  as  anomalous  earthquakes.  In 
addition,  there  was  one  event.  No.  43,  for  which  a  Love-wave  spectrum  could 
not  be  determined,  which  was  anomalous. 
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The  following  conclusions  ere  noted  from  this  study: 

Discrimination  between  earthquakes  and  underground  nuclear  explosions 
in  Eurasia  by  the  m^ :Mg  criterion  is  enhanced  when  the  amplitude  of  Love 
waves  at  a  period  near  20-sec,  rather  than  Rayleigh  waves  (vertical- 
component),  is  used  to  determine  Mc.  In  general,  Love-wave  magnitude  can 
be  used  to  reduce  the  number  of  anomalous  events.  Of  nine  Eurasian  events 
with  anomalous  m^:M^  values,  six  are  also  anomalous  earthquakes  by  the 
m^:M^  criterion.  None  of  the  14  Kashmir  mainshock-aftershock  sequence 
earthquakes  has  been  found  to  be  anomalous  although  4  events  in  the 
sequence  have  anomalously  low  Rayleigh-wave  magnitudes. 

The  Love-wave  amplitude  spectrum  (A^)  also  appears  to  be  a  useful 
earthquake-explosion  discriminant.  By  the  m^ :A^  criterion,  only  four  of 
the  Eurasian  earthquakes  studied  were  anomalous.  For  the  remaining  two 
events  (events  No.  37  and  43)  with  anomalous  m^M^  values,  the  shape  of 
the  Love-wave  amplitude  spectrum  of  event  No.  37  has  normal  shallow  earth¬ 
quake  characteristics.  There  are  no  available  spectral  data  for  event 
No.  43,  because  the  surface-wave  motion  was  too  small. 

For  a  Kashmir  mainshock-aftershock  sequence  studied,  the  focal 
depths  and  focal  mechanisms  were  the  same  for  almost  all  events.  Furthermore 
by  both  the  and  :Ag  criteria,  all  the  events  in  the  sequence  are 

non-anomalous.  In  summary,  it  is  concluded  that  Love  waves  provide  a 
powerful  earthquake-explosion  discriminant. 

Although  others  have  suggested  that  anomalous  earthquakes  are 
confined  to  the  interior  of  the  continent  or  to  aftershock  sequences,  two 
oceanic  events  and  none  of  the  Kashmir  aftershock  sequence  are  found  to  be 
anomalous  earthquakes  by  the  Love-wave  discriminant  employed  in  this 
study.  It  is  believed  that  most  anomalous  earthquakes  can  be  reclassified 


as  non-anomalous  by  an  adequate  analysis  of  the  available  data. 
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ANOMALOUS  EARTHQUAKES  STUDIED 
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The  speetra  were  all  normalized  to  an  epieentral  distance  of  9 


Figure  l.  versus  values  for  nine  Eurasian  events  studied.  The  number  given 
in  parentheses  is  the  number  of  observations. 
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Figure  **•  Fault-plane  solution  for  event  KL  of  the 
nor  thv;e  stern  Kashmir  nainshock-af ter  shock 
sequence*  The  octagon  and  triangles 
stand  for  compressional  and  dilatational 
P-wave  motion,  respectively* 
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Figure  5«.  Amplitude  spectrum  of  fundamental-mode  Love 
waves  for  event  No*  37  reoorded  at  VLPE 
station  KON-,  The  peak  magnification  for 
the  high  gain  instrument  (solid  curve)  is 
about  20  times  that  for  the  low  gain 
instrument  (dashed  curve)* 
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Figure  5d»  Amplitude  spectrum  of  fundamental-mode 
Love  waves  for  event  No.  135* 
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Figure  5e.  Amplitude  spectra  of  fundamental-mode 
Lore  waves  for  event  No.  145* 
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Figure  5f.  Amplitude  speotra  of  fundamental -mode 
Love  waves  for  event  No*  146. 
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Figure  6a.  Amplitude  spectra  of  fundamental -mode 

Love  vraves  for  Kashmir  earthquakes  record 
ed  at  SHI  (Shiraz,  Iran).  The  numbers 
following  the  station  code  are  epicentral 
distance  and  azimuth  to  the  station.  The 
character  attached  to  each  ctirve  repre¬ 
sents  the  particular  Kashmir  event. 


AMP  (CM-SEC) 


QUE  (7.9°,  225.6°) 


10  "N 


-5 


1  0  -J — - — i i  «  i— r-r-rrj — - — i 1 — r-nrrrr 

10°  10'  10 

PERIOD  (SEC) 


Figure  6c.  Amplitude  spectra  of  fundamental -mode 

Love  waves  for  Kashmir  earthquakes  record 
ed  at  QUE  (Quetta,  Pakistan). 
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figure  6d.  Amplitude  spectra  of  fundamental -mode 

Love  waves  for  Kashmir  earthquakes  record' 


▼alues 


MODELLING  COMPLEX  EARTHQUAKES 
by 

Jaime  Yamamoto 

Introduction 

It  la  well  known  among  seismologists  that  seismograms  from  large 
shallow  earthquakes  very  often  have  a  complicated  appearance  near  the 
beginning  of  the  record.  This  complex  appearance  presents  itself  as  small 
amplitude  onsets  prior  to  the  P-wave  arrival  and  sometimes  as  small 
glitches  on  the  P-wave  trace  which  noticeably  distort,  its  shape.  The 
main  characteristic  of  these  onsets  and  glitches  is  that  in  general  they 
cannot  be  associated  with  any  known  phase. 

The  complication  in  shallow  focus  earthquake  records  may  be  explained 
in  at  least  four  ways:  (1)  structure  beneath  the  epicenter  and/or 
interactions  with  the  free  surface,  (2)  earthquake  multiplicity,  (3)  ray- 
path  effects  and  (4)  station  effects. 

e 

The  last  two  effects  mentioned  in  the  above  paragraph  can  be  diminished 
by  studying  records  from  stations  at  an  appropriate  range  of  epicentral 
distances  (30a  -  80*)  and  by  observing  only  features  that  can  be  observed 
at  several  stations.  Therefore,  our  main  concern  at  this  moment  is  to 
investigate  whether  the  complication  shown  by  some  shallow  focus  earthquakes 
is  produced  by  interaction  with  the  structure  in  which  the  source  is 
embedded  or  is  a  multiplicity  effect. 

Method  of  Analysis 

The  technique  will  be  to  match  the  observed  long-period  teleseisralc 
P-wave  signals  with  synthetic  seismograms  computed  by  the  Uaskell-Thomson 


73 


method.  The  program  HASKELL  written  by  Herrmann  (unpublished,  1974)  with 
slight  modifications  was  used  for  the  generation  of  the  required  synthetic 
seismograms.  The  theory  used  in  writing  this  program  is  contained  in 
Haskell  (1964)  and  Hudson  (1969)  and  an  application  appears  in  Herrmann  (1976) 

In  this  program  the  far-field  source  time  function  is  represented  by 
a  trapezoid  of  area  numerically  equal  to  the  seismic  moment  of  the 
earthquake.  The  range  of  epicentral  distances  is  restricted  to  30-80*  to 
avoid  the  complications  produced  by  the  heterogeneities  of  the  upper-most 
mantle.  The  program  options  include  source  crust  response,  receiver  crust 
response,  anelastic  attenuation,  and  Instrument  response. 

Some  Results 

The  method  is  being  applied  to  several  earthquakes  which  have  occurred 
near  the  coast  of  Chiapas,  Mexico.  In  this  progress  report  I  will  present 
some  preliminary  results  for  two  earthquakes  that  occurred  in  that  area 
on  April  29,  1970  at  11:22  GMT  (Mg  -  6.3)  and  April  30  at  08:33  (Mg  -  6.4). 
The  epicentral  coordinates  are  14.50N,  92.77W  and  14.58H,  93.20W,  for  the 

e 

first  and  second  event  respectively.  Both  earthquakes  have  almost  the 
same  depth  of  about  15  km  based  on  pP  -  P  readings.  From  P-wave  first  motion 
data  (54  compressions  and  1  dilatation) ,  the  two  earthquakes  show  identical 
fault-plane  solutions  (strike  ■  N92E,  dip  ■  20*,  slip  -  71°).  The  mechanism 
represents  reverse  faulting.  The  spatial  distribution  of  aftershocks  was 
used  as  an  aid  in  deciding  which  of  two  nodal  planes  correspond  to  the 
fault  plane.  It  should  be  noticed  that  the  nodal  planes  are  not  well 
constrained  due  to  the  fact  that  only  one  station  recorded  a  dilatation. 

In  future  studies,  I  plan  to  use  other  methods  in  an  attempt  to  improve 
the  solution. 


The  two  events  are  part  of  an  earthquake  sequence  that  took  place 
about  100  km  off  the  coast  of  Chiapas,  Mexico  in  April,  1970.  They  were 
chosen  because  they  show  similarities  in  magnitude,  location  and  source 
geometry,  but  at  the  same  time  the  appearance  of  the  long-period  P-wave 
signals  are  quite  different  (see  Figures  1  and  2). 

The  shape  of  the  P-wave  for  the  earthquake  which  occurred  at  08:33 
as  observed  at  the  station  BK.S  (distance  ■  35  degrees,  azimuth  «  317  degrees) 

looks  very  simple  suggesting  a  simple  rupture  process.  The  synthetic  seis¬ 
mogram  displayed  on  the  same  figure  was  computed  assuming  es  a  source 

function  an  isosceles  triangle  with  a  base  of  5  seconds  and  an  area 

25 

arbitrarily  taken  equal  to  10  .  The  velocity  structure  at  the  source 

region  was  obtained  from  extensive  refraction  studies  carried  out  in  that 
area  (Shor  and  Fisher,  1961).  The  structure  beneath  the  recording  station 
was  assumed  to  be  a  simple  two-layer  crust  typical  of  continental  regions. 

The  agreement  between  the  synthetic  and  observed  seismograms  in 
Figure  1  is  quite  good  for  the  first  30  seconds.  For  later  times  the 
seismogram  is  largely  controlled  by  reflections  and  conversions  coming 
from  deep  Interfaces.  At  these  depths  (about  200  km)  the  velocity  model 
is  not  known  very  well. 

The  observed  record  of  the  11:22  earthquake  exhibits  a  far  more 
complicated  character  (see  Figure  2) .  It  shows  an  extra  secondary  arrival 
superimposed  on  the  first  down  swing.  The  amplitude  of  the  second  upswing 
has  been  severely  attenuated.  The  relative  amplitude  of  a  third  upswing 
has,  on  the  other  hand,  increased.  The  features  are  common  to  most  of 
the  stations. 

A  set  of  synthetic  seismograms  were  generated  assuming  a  simple 
source  time  function  located  at  various  depths.  None  of  these  synthetics 
look  like  the  observed  seismogram.  Therefore,  we  have  to  resort  to 


multiple  events  in  order  to  explain  the  observed  complexities.  A  source 
function  composed  of  two  isosceles  triangles  of  base  equal  to  5  seconds 
and  separation  between  peaks,  equal  to  5  seconds  was  used  to  simulate  two 
consecutive  events.  The  seismic  moment  of  the  second  event  was  chosen  to 
be  two  times  larger.  The  synthetic  and  observed  seismograms  are  shown 
on  Figure  2.  The  first  13  seconds  of  seismogram  agree  reasonably  well. 

An  extra  second  arrival  approximately  5  seconds  after  the  first  P  arrival 
is  clearly  shown  by  the  synthetic.  For  greater  times,  there  is  a 
considerable  discrepancy  between  the  theoretical  and  observed  signals. 

In  particular,  a  large  second  upswing  is  predicted,  but  is  not  seen  on 
the  real  seismogram.  The  computed  third  upswing,  on  the  other  hand,  is 
small  as  compared  with  that  of  the  observed  record.  At  this  time,  the 
reasons  for  the  discrepancy  are  not  known.  The  rupture  process  of  this 
earthquake  is  likely  to  be  more  complicated  than  assumed,  therefore  the 
simple  superposition  of  two  sources  may  not  work  properly.  Further 
work  will  concentrate  on  using  more  complex  source  models. 
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Figure  1.  Synthetic  (left)  and  observed  (right)  P-wave  signals 
at  BKS  (distance  ■  35*  and  azimuth  *  317*)  for  the 
April  30,  1970  event.  Synthetic  was  computed  for  a 
source  depth  of  15  km.  The  amplitude  units  are  arbitrary 
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Figure  2.  Synthetic  (left)  and  observed  (right)  P-wave  signals  at 
BKS  for  the  April  29,  1970  event.  Synthetic  computed 
for  a  source  depth  of  15  km.  .  The  amplitude  units  are 
arbitrary. 
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THE  EFFECT  OF  COMPLEX  NEAR-SOURCE  CRUSTAL  STRUCTURE 
ON  THE  WAVE  FORMS  OF  RAYLEIGH  WAVES 


Brian  J.  Mitchell  and  Chiung-Chuan  Chang 


An  extensive  study  of  the  attenuation  of  Rayleigh  wave  spectral 
amplitudes  across  Eurasia  has  been  reported  by  Yacoub  and  Mitchell  (1977) . 

The  confidence  limits  of  their  observed  Rayleigh  wave  attenuation  coefficients 
were  substantially  larger  than  confidence  limits  for  attenuation  coefficients 
in  North  America  (Herrmann  and  Mitchell,  1975).  The  large  confidence 
limits  for  the  Eurasian  data  probably  reflect  the  complex  structure  across 
which  most  of  the  Rayleigh  waves  passed.  Lateral  refraction  and  multi- 
pa thing  which  originate  along  the  propagation  paths  of  the  Rayleigh  waves 
are  obvious  possibilities  as  causes  of  scatter  in  the  data.  Another 
possibility  is  that  lateral  complexities  in  the  source  region  perturb  the 
Rayleigh  waves  as  they  are  being  formed.  The  purpose  of  the  present  research 
is  to  investigate  the  extent  to  which  near-source  complexities  affect  the 
wave  forms  of  fundamental  mode  Rayleigh  waves  observed  at  distant  stations 


in  Eurasia. 


Methodology 


The  effects  of  near-source  structure  and  nulti-pathing  on  the  wave 
forms  of  Rayleigh  waves  can  be  studied  by  comparing  the  observed  wave  forms 
with  synthetic  Rayleigh  wave  seismograms.  Dr.  R.3.  Herrmann  has  written 
a  program  for  such  computations  which  he  has  permitted  us  to  use  for  this 
study.  Necessary  input  for  theoretical  seismogram  calculations  are  the 
following:  (1)  the  fault-plane  solution  for  the  earthquake  which  produced 
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the  Rayleigh  waves  being  studied,  (2)  the  amplitude  spectrum  at  the  source, 

(3)  attenuation  coefficient  values  as  a  function  of  period  for  each  mode, 
and  (4)  a  knowledge  of  the  instrument  response  for  the  seismographs  used. 

Fault-plane  solutions  for  these  computations  are  available  for  three 
of  the  earthquakes  studied  by  Yacoub  and  Mitchell  (1977) .  Two  of  these 
have  been  used  as  sources  for  synthetic  seismogram  computations  in  this 
study.  In  addition,  synthetic  seismograms  were  computed  for  a  nuclear  event 
at  Novaya  Zemlya.  The  source  spectrum  was  estimated  from  the  magnitudes 
of  the  events  and  the  scaling  relationships  of  Aki  (1967).  Attenuation 
coefficients  corresponding  to  shield  regions  for  fundamental  and  higher 
modes  have  been  computed  by  Mitchell  et  al  (1977)  and  were  used  for  all 
computations  in  this  study.  Instrument  responses  were  computed  for  the 
long-period  WKSSN  instruments  used  for  this  research.  Seven  modes  were 
used  for  all  of  the  computations  in  this  work.  These  were  computed  for 
the  Baltic  shield  model  of  Noponen  et  al  (1967) . 

Results 

# 

The  synthetic  seismograms  calculated  as  described  above  were 
compared  to  observed  seismograms  for  selected  events  and  stations.  Figure  1 
(from  Yacoub  and  Mitchell,  1977)  is  a  map  of  the  events  and  stations  for 
which  observed  Rayleigh  wave  seismograms  are  available.  Synthetic  seismograms 
were  first  computed  for  events  near  a  transition  between  stable  and  tectonic 
regions,  and  for  paths  which  were  entirely  across  stable  regions.  The 
earthquakes  with  known  fault-plane  solutions  which  occurred  near  such  a 
transition  arc  those  of  October  28,  1971  and  November  18,  1971.  Stations 
to  which  paths  were  restricted  to  stable  regions  are  those  in  northern  Europe. 
Results  of  the  first  computations  appear  in  Figures  2  and  3.  The  agreement 
between  the  main  features  of  the  observed  and  theoretical  seismograms  is 


is  good.  The  durations  and  shapes  of  the  envelope  of  the  observed  wave 
forms  are  similar  to  those  for  the  computed  wave  forms  in  all  cases 
except  that  of  the  November  18  event  as  recorded  at  NUR.  The  duration 
of  the  envelope  of  the  first  100  seconds  is,  however,  very  similar  to  that 
of  the  synthetic  seismogram.  Later  energy  bursts  appear  on  the  observed 
seismogram  which  are  not  predicted  theoretically.  These  cannot  be  caused 
by  near-source  complexities  since  they  are  so  much  later  than  the  first 
energy  group;  they  may  be  due  to  multi-bathing  which  arises  along  the 
propagation  paths.  Although  near-source  complexities,  for  the  above 
examples,  do  not  have  a  large  effect  on  the  gross  features  of  the  wave 
forms  studied,  they  may  have  a  perturbing  effect  which  is  difficult  to 
separate  from  other  effects,  such  as  multl-pathing  or  departures  of  fault- 
motion  from  that  of  Che  idealized  sources  which  are  assumed  for  the 
synthetic  seismogram  computations. 

A  Novaya  Zemlya  nuclear  event  (November  7,  1968)  was  studied  as  an 
example  of  a  source  near  a  continental  boundary.  Synthetic  seismograms 
'were  computed  in  the  same  way  that  they  were  for  earthquakes,  except 
that  the  source  radiation  pattern  was  taken  to  be  circular  and  the  source 
spectrum  Included  a  peak  such  as  that  which  appears  in  published 
source  spectra  of  nuclear  explosions  (e.g.  von  Seggern  and  Lambert  (1970)) 

Observed  and  theoretical  Rayleigh  waves  from  the  Novaya  Zemlya 
event  of  November  7,  196C  to  stations  MSI1  and  KliL  appear  in  Figure  4 

The  observed  seismograms  differ  substantially  from  the  theoretically  pre¬ 
dicted  seismograms  in  all  cases.  The  reason  for  the  differences  cannot 
as  yet  be  ascertained.  Possible  reasons  for  the  differences  include  the 
effect  .teral  complexities  near  the  source,  the  presence  of  a  water 
layer  near  the  source,  or  departures  of  the  source  spectrum  from  that 


which  was  assumed 
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Rayleigh  waves  for  paths  across  the  tectonic  regions  of  Figure  1 
differ  greatly  from  those  which  cross  the  stable  regions  (see  Figure  5) . 
Different  sub-regions  within  the  broad  region  designated  as  tectonic  also 
produce  seismograms  which  differ  considerably  from  one  another.  For 
instance,  the  seismogram  recorded  at  a  station  in  the  Mediterranean  region 
(AQU)  has  a  very  poorly  developed  Rayleigh  wave  which  appears  to  be 
dominated  by  multi-pathing.  The  other  two  stations  (MSH  and  SUL)  in 
central  Asia  have  well-dispersed  wave  trains  which  differ  greatly  form  those 
recorded  for  shield  paths.  Further  synthetic  seismogram  computations 
are  in  progress  for  models  corresponding  to  the  tectonic  regions  of 
Eurasia.  These  should  tell  us  if  the  features  of  the  observed  seismograms 
are  produced  by  the  velocity  structure  in  the  region,  or  if  multi-pathing 
or  near-source  complexity  are  important. 
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Figure  2.  Observed  and  synthetic  Rayleigh  wave  seismograms  for  the  USSR 
earthquake  of  October  28,  1971.  The  Rayleigh  waves  travelled 
almost  entirely  across  stable  regions  .0  stations  KEV  and  NUR 
The  time  scales  Indicate  the  duration  of  the  plotted  trace  in 
minutes . 
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Figure  3.  Observed  and  synthetic  Rayleigh  wave  seismograms  for  the  USSR 
earthquake  of  November  18,  1971.  The  Rayleigh  waves  travelled 
almost  entirely  across  stable  regions  to  stations  KEV  and  NUR. 
The  time  scales  indicate  the  duration  of  the  plotted  trace  in 
minutes. 
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Figure  A.  Observed  and  synthetic  Rayleigh  wave  seismograms  for  the  Novaya 
Zemlya  nuclear  event  of  November  7,  1968.  The  Rayleigh  waves 
travelled  almost  entirely  across  stable  regions  to  stations 
MSU  and  KbL.  The  time  scales  Indicate  the  duration  of  the  plotted 
trace  in  minutes. 
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Figure  5.  Observed  seismograms  from  the  USSR  earthquake  of  November  18,  1971 
The  Rayleigh  waves  travelled  along  complex  tectonic  paths  to  the 
stations  AQU,  MSH,  and  SHL.  The  time  scales  Indicate  the  duration 
of  the  plotted  trace  In  minutes. 
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